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a^ : 

0\ . For the quark-gluon plasma, an energy-momentum tensor is found corresponding to 

^ \ the high-temperature Braaten-Pisarski effective action. The tensor is found by considering 

the interaction of the plasma with a weak gravitational field and the positivity of the 

■ energy is studied. In addition, the complete effective action in curved spacetime is written 

. down. 
0\ ■ 

(N ; 
T— I ! 

>■ . 1 Introduction 

^ ! The high-temperature properties of the quark-gluon plasma are of some interest, both 

^ \ in their own right and as a starting point of a re-summation of perturbation theory [1]. 

^ \ These properties are encoded in an 'effective action' F, which has a coefficient proportional 

^ ■ to (where T is the temperature) [2]. This effective action is nonlocal, but in a very 

G^ ■ special way. 

^ . If F is taken rather literally as an action for an effective field theory, a question 

I \ which comes to mind is: what is the associated energy-momentum tensor density T^jj[x) 
and energy-momentum vector P^l Since F is nonlocal (in time as well as in space), the 

' canonical construction for T^jj is not available. An alternative definition of T^jj is given 

>• ■ as follows. First generalize F to a curved spacetime background (which is asymptotically 

* I — 

^ , Minkowskian) , giving a functional F of the vierbein (e^") and metric (as well as of the 
^ \ gluon and quark fields). Then 



T^a (•^) 
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In order to find T^a in this way, we need to calculate hard thermal loops with one graviton 
line emerging in addition to the quark and gluon lines. In Section 2 we perform this 
calculation for the quark-antiquark- graviton case and for the two-gluon-graviton case. 

The calculation of hard thermal loops always gives results which may be written in 
the form 

9'^7:rTs / dnx{Q^), (2) 
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where is a lightlike 4-vector and dfl is an integral over the directions of Q. Then X 
is not unique, but different choices for X have different properties. The best choice may 
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depend upon the use to which F is being put. In this paper, we make the choice in which 
X is Lorentz invariant (though of course (2) is not) [3]. Our experience with hard thermal 
loops indicates that this X enjoys the following properties:- 

(a) It is nonlocal, but all the nonlocalities come from products of {Q.d)~^ operators. 

(b) It is homogeneous in Q of degree zero, and has dimensions {energy)~'^. 

(c) It is (as already stated) Lorentz invariant, with respect to Lorentz transformations 
in the asymptotic Minkowski space. 

(d) It is explicitly gauge-invariant (not just BRST invariant). 

(e) In the case of f, it is invariant under general coordinate transformations, Weyl 
transformations, and (when spinors and vierbeins occur) local Lorentz transformations, all 
of which are restricted to tend to the identity at infinity. 

We conjecture that these properties are sufficient to to fix T uniquely, given only the 
lowest order terms. In Appendix B we give an example of this. In Section 3, we write down 
a all-orders expression for f which satisfies the above 5 conditions and which, therefore, 
we believe is the unique correct result. To one graviton order, it agrees with the explicit 
calculation in Section 2. 

Each of the properties in (e) gives rise to a Ward identity which T^^ should satisfy, 
and which we check for the tensor calculated in Section 2. 

Since f is nonlocal, T^a is also; so it is not really a density. However the energy- 
momentum is given by 



One of the Ward identities ensures that, provided we use the equations of motion given by 
varying F, Pa is independent of time. (T^^ is also symmetric and traceless if the equations 
of motion are used.) 

The obvious questions to ask about Pq concern positivity. We have found no positivity 
property for Pq in general, but we do find that it is positive if the equation of motion (to 
one- loop order) is used. But this latter property tests very little about Too: it actually just 
tests the positivity of the residue of the propagator at the pole (whose position is given by 
the solution of the equation of motion). This is something that could have been checked 
without knowing T^^,. 

In [4] , we have previously discussed all these points for the model of scalar field theory. 

2. The calculation of T^i^ from hard thermal loops. 

In this section we use thermal field theory to calculate X in (2) to zeroth order in the 
QCD coupling constant g. 

Our notation is as follows. We take QCD with an SU{N) colour group and Nf 
flavours. The quark field is -0, the gluon field with A,B, .. — 1, (N"^ — 1). The 
vierbein is e^a, with i^, ... being Riemann indices and a,b, ... being local Lorentz indices, 
are hermitean colour matrices. We define 



The graphs are in Fig 1, and the Feynman rules are listed in Appendix A. The calcu- 
lation of hard thermal loop is fairly standard, and we will just state the result. (We used 
Mathematica to do the algebra.) 




(3) 



Cg = N, Nftx{tH^) = CgJ^-^, t^"^ = C'qI. 



(4) 
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We write 

[{Cg + C,)A^^{y')V^:{T,p)At{y) + C'^i^{y')SUr.p)m] ■ (5) 

Then 

V^; = QaQpQ^Q'' {P^K.R + \K^R.P) - QaQ^'Q" {^R^Pp + R.PK^) 
-QaQpQ" {K.RP" + K^R" + R'^K") + Q^Qp [R'^K^ + i^^^^.^ ^ Q^Q^R^Pp 
+2Q^Q^ {RpP^ + R^Kp + R%) - (2P''rj^ + K^r/'^^) - 2Q^i?^< + <r/J 

+ (p, z/) ^ (r, //) + (a^P), (6) 

where 

Hip + Oa^M + + ^)f^Qa - (P + r).KQaQ^Q ■ 7} 

+ ^A;.Q{a,^, g.7} + iQJfcVb^, Q.7} + iQ^I^'^ba' Q-T}] (8) 

(where a ah = i(7a76 - 7&7a))- 

There are certain Ward identities which express, to lowest order, the invariance prop- 
erties of r. Invariance under local coordinate transformations implies 

2A:^^aT = Si;:kxU>^''{p)+paU^^''{r) - {{p,u) ^ (r,/.)), (9) 

k^Sxa^ra^ip)-Pa^ir), (10) 

where 11 and S are the gluon and quark 2-point functions from hard thermal loops [1] 

S(p) = 3§^, n'^^ = 4 (r)"'' - P^Q"" - P'^Q'^ + P^Q'^Q'') (11) 
Q.p 

(in the notation (7)). 

Invariance under local Lorentz transformations implies the identity 

Sxa - SaX = (TXa^ip) - S(r)(7Aa. (12) 

Finally, QCD gauge invariance gives 

P^^Vr, = r,V^; = 0- (13) 
3 



In hard thermal loops, quark masses are neglected, and we expect Weyl invariance, 
since this holds for massless QCD. This is invariance under 

Ai ^ Ai, e^a - e-^^^^e^a. g^u - e-''<^^g^., ^ e^^^^^. (14) 

It implies the identities 

ry^'^Vr = 0, (15) 
r7^«5Aa = i[S(p) + S(r)]. (16) 
All the above identities are satisfied by (6) and (8). 

3. Hot QCD in curved space 

In this section, we guess a form for the effective action, f , in curved space, to all orders 
in g and all orders in the gravitational field. Our method is simply to guess a functional 
of the quark field the gluon field A, and the metric and vierbein, which satisfies all five 
conditions (a) to (e) of Section 1. We believe that this functional is likely to be unique, 
because of the sort of arguments used in Appendix B. 

To write down f we must introduce some notation. Let y^{x, 9) be the null geodesic 
satisfying 

y{x, 0) = a;, ^ ^ Q as ^ ^ -oo. (17) 
As in [5], we assume that spacetime is asymptotically Minkowskian, and that 

^ ^ g as ^ ^ +00 (18) 
also. The status of this assumption is explained in [2,5]. Define 

[/(x;^2,^i) = Pexp^ r de'y.A[y{x,e% (19) 

where A and U are here matrices in the adjoint representation of SU (N), A being defined 
by 

^BC ^ fABCj^A_ (20) 

P denotes path-ordering with respect to this matrix multiplication and dot denotes differ- 
entiation with respect to d' . 
Define A by the equation 

A{x)= / d^e^[2'(^'^)] / d9'e-^^y^''^^'^^{R^Mx,0'Wy'' + Ay2}, (21) 

Jo J-oo 

where here R is the Riemann tensor. This equation defines A implicitly: for example it 
may be used to generate a series for A in powers of the gravitational coupling n (since R 



is of order k). The importance of A is that, under the Weyl transform (14), it transforms 
as 

SA{x) = 2a{x). (22) 

The proof of this uses 

= -2[a + &% (23) 

and the fact that 

^0^A[y(x,e)] (24) 

is Weyl invariant. 

With these definitions, the gluon part of f is of the form (2) with 

Xg = ^^26^^'^ / ^'^A/=^^'''(^)e^^'^tr[L»^t/(a;;oo,0)i;,t/(a;;0,-oo)], (25) 

where D is the colour covariant derivative 

(D^C/)^^ = d^U^"" - gA^'^U^''. (26) 

(The factor Co in the denominator in (25) is due to our choice to write the colour matrices 
in the adjoint representation.) To check that (25) satisfies the 5 conditions in Section 
1, we note that the expansions for y in terms of k and for U in terms of g and k, create 
nonlocalities of the required kind, which are homogeneous of degree zero in Q (see examples 
in [3] and [5]). Invariance under general coordinate transformations is obvious since U is 
invariant under such transformations. The factor fixes up Weyl invariance. 
To leading order in g we have 



D^JJ{x\^, -oo) ~ g 



J —Oi 



dO y.A 



(27) 



(28) 



and hence one may verify that (25) gives (6) to this order and to first order in n. 
In fiat space 

D^U{x; 0, -oo) = -9 d9 U{x; 0, 9)Q^Fx^{x + Q9)U{x; 9, -oo), 

J —oo 

and hence (25) reduces to 

_i_ rO nOO 

^r^g" \ d9 / d9\Y{Q^Fx^,{x + Q9')U{x; 9', 9)Q,F''^'{x + Q9)U{x; 9, 9')] 
20g Jo 

/o j-e 
d9 d9'Q^F^^{x + Q9')[U{x;9',9)]'^''Q,F''^''{x + Q9). (29) 
-oo J — oo 
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where we have used the same type of arguments as in the first reference [2]. The second 
form in (29) is in agreement with [3]. 

For the quark part of the effective action, we need one further definition. Let Qx be 
the 'spin connection' 

n, = -^e'^e%,aab. (30) 

Also let (with A defined in (21)) 

= Q,- ie«e^^A,^aa6, (31) 

which is a Wcyl invariant, but retains the same property as Q under local Lorentz trans- 
formations. Then define 

W{x;e2,9i)=Pexp [ ' den'x[yix,e)]y^. (32) 

This has the same relation to local Lorentz spinor transformations that U in (19) has to 
local SU (N) colour transformations. 
Define also 

*(a;)=/ dee^^^y^'^'^^^U{x;0,e)W{x;0,e)i;[y{x,e)], (33) 

^— oo 

where U is now in the fundamental representation of SU{N) and acts on the colour index 
of i/j just as W acts on the Dirac index. Because of the transformation property of i/j in 
(14) and of (24), \& is Weyl invariant as well as transforming the same way as ip under 
local colour and Lorentz transformations. 

With these definitions, the quark part of f is given by (2) with 

Xq = \iC'q J d^x^y-g{x)ei^^''^[y%=oea^{x)i;{x)'y''<^/{x) + herm. conj.. (34) 

Here the e^^ fixes up Weyl invariance (see (14) and (22)). 
One may verify that (34) gives (8) to zeroth order in g. 

We are struck by the fact that these effective actions, although they are nonlocal, seem 
to be just as unique as the ordinary QCD action in curved spacetime. This is because the 
extra freedom allowed by the special kind of nonlocality (property (a) in Section 1) is 
compensated for by the different dimensionality (property (b)). (See Appendix B for an 
example.) The invariance properties (c), (d) and (e) are common to the two cases. 

4. Positivity properties 

The energy-momentum tensor of an ordinary local classical field has an energy-density 
Too which is a local positive definite functional of the fields. What, if any, are the cor- 
responding properties of (5), (6) and (8)? Since Too is not local there can be no energy- 
density; so we study the total energy Pq given by (3). Then we require k = in (6), (7) 
and (8), and K^o in (7) reduces to S^q. 



6 



Let us test the positivity of the gluon contribution to (3) by inserting into (5) the 
particular field 

Ai = e^cos{s.y), (35) 

where s is an arbitrary 4-vector and e satisfies Cq = s.e-^ = 0. The contribution to (3) 
thus obtained is proportional to 



efeMT'S^O), (36) 



which is proportional to 



£6 
s2 



.0 iJfo±hl_2 
So - f 



(37) 



This expression is unambiguous for sq > |s|. Like most expressions in thermal field, it has 
a branch cut for sq < |s|, and requires further definition there. But it is clear that, to test 
the positivity of the energy, we should take the real part. It is not positive for all values 
of so/|s|; so we have a counter-example showing that Pq is not a positive functional of the 
fields. 

The next thing we can do is study the effect of imposing the restriction that the 
fields in (5) should satisfy the equations of motion given by the free action plus F to one- 
loop order. (In the Braaten-Pisarski resummation method the one-loop g^T^ term in (5) is 
considered to be effectively of the same order as the free action.) Then Pq is independent 
of time in virtue of the Ward identities (9) and (10). 

But if we do impose the equations of motion it turns out that these Ward identities 
are sufficient to determine Pq: we do not actually need (6) and (8). Take for example the 
gluon self-energy. Define the coefficient 

c=-^{Cg + C,), (38) 

and the transverse and longitudinal tensors (actually the negatives of projection operators 
in our metric with tjqq = +1) 

= i^ij - PiPj)maVjl3, Paf3 = 'VafS + - Paf3 ih j = 1, 2, 3). (39) 

We call the complete self-energy function (free plus thermal) 11, and write 



n«;3 = n^pj. + n^pi'., (40) 



where 

= -cn^ + p^ = -cn^ + (4i) 

Similarly we denote the complete energy function by Vqq . There is a Ward identity 
connecting V with II, completely analogous to (9). In this identity, set k = and the 
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suffix a = 0, differentiate with respect to ko for fixed po, and insert (40) into the righthand 
side, to obtain 



2 = P^^"' 



Now use the field equations 



Po 



d 

dpo 



Po 



d 2pl 



dpo p2 



+ 1 



= 0, = 0. 



(42) 



(43) 



Then the question is about the positivity of the differential functions on the right of (42) 
at the solutions of (43). These are nothing but the residues of the poles of the propagator 
given by the solutions of (43). 

We will treat the gluon energy case as an example. Define 



Po 



Then [1,2] 
So 



x=P^, L = x\n\{x + l)/{x-l)\. 



-p2(l - x^) - c[2x^ + (1 - x'^)L] = 0. 



d - d " 

2 V^o = x—U^ = x—U^ - 2U^ = 2p^ + c[(l + x'^)L - 2x^\ > 0. 



For the longitudinal part, we have 

= (l-a;2)[p2 + 2c(2-L)], 

and we find that, at fl^ = 

Foo = c[2a;2 + (1 - x'^)L\ > 0. 



(44) 

(45) 
(46) 

(47) 
(48) 



The quark energy may be treated similarly and the result is also positive. 

Thus the only positivity property for Pq which we have been able to find is nothing 
more than the positivity of the residue at the pole of the efi^ective propagator. 

Appendix A 

In this appendix we present the QCD Feynman rules in a curved spacetime. The 
Lagrangian is: 

^ = ^yM + ^fix + ^ghost + I^Ti (^1) 



where 



< = Q^^t - d.A^ + ^/^^^AjA^; [A = 1, 2, • • • , AT^ - 1] , (^2) 
Cf,, = -\ ( VM;^) i^^At), (A3) 



g^^- = r]^- + K^-- ^-det{g) = 1 - '^rj^.r" + 0{k''); [k = ^32^] (^6) 
V^A^ = d^A^ + r^^^«; = \g^P {d^gp^ + d^gp^ - dpg^^) {Al) 



= + ^^a6<, (^9) 
<^ab = ^ (7a 76 - 767a) , (-410) 

< = ^e"'^ (a^et - a.ey + le^^e^'^ (a.e^ - a,e^) - (a ^ h). {All) 

The equations above show that the graviton fields (j)^^ couple to the gluon fields A^ and the 
ghost fields The vierbein h^i, couple to the quark fields tfj and there is also a coupling 
with the gluon field and the quark field. The form of these interactions is such that the 
theory is invariant under BRST transformations, local coordinate transformations, local 
Lorentz transformations and the Weyl transformations given by Eq. (14). 

From the Lagrangian in Eq. (Al), one can now obtain the momentum space gauge 
particles propagators and the couplings with one graviton. In all the expressions which 
follows we will always denote the graviton indices by (ct, (3) (or (a, h) when quarks are 
involved). We reserve fi, v, p and a for the gluons. 

The Feynman rules are: 

gluon propagator 

5^^ ^ {A12) 
two gluons one graviton coupling 

nVctjS {kfx -\- kj^ + 2kn k^, — k • k T]fj,ir) 



k^Vfj.f ^^VtSu kj^ k^rjoifx 2,k^ k^ rj^j^ + k • k^r}oi^r}i3,j 
+ symmetrization under (a ^ /?)] 

+ {k'^,li,A) ^ {k^,v,B) {AlZ) 
three gluons coupling 



three gluons one graviton coupling 



2 

+permutations of (A;^,//,^), {k^,v,B), (k'^,p,C) {A15) 
four gluon coupling 

[(/^^'^^ - /^^'^^) V,.Vp. + (/^^'^^ - /^^'^^) V,pV.. 

jAB,CD _ jGAB jGCD {AlQ) 

four gluon one graviton coupling 

2'^ nAB,CD fVupVi^aVa/d \ 

+permutations of (/c^,//, A), {k^,u,B), {k^,p,C), {k^,a,D) {A17) 

ghost propagator 

ghost-ghost graviton coupling 

^S^^ {r^^.k^ . k^ - k^k^, - k^k^) iA19) 

ghost-ghost gluon coupling 

igf^^^k^ (A20) 
ghost-ghost gluon graviton coupling 

/^''^ ik^V(3, + k^Va, - Vapk^) {A21) 

quark propagator 
-5.,^; [i, j = l,2,--.,N] {A22) 
quark-quark vierbein coupling 

^ij ^ [2 {P' - P)a lb - '^riabl ' {p' - p) - 1 ' {p + p') (Tab - (^abl ' {P + p')] {A23) 
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quark-quark gluon coupling 

91^.T,^ {A24) 
quark-quark gluon vierbein coupling 

9 « iVai^lb - Vabli^) T^j (^25) 

In the expressions above, each gluon is labeled by (A;"^, A)^ {k^ , u, B), {k^ , p, C) 
and {k^ , cr, D), where k^ is the gluon momentum. The ghosts label are {k^ , B) and 
C), where k^ is the ghost momentum. The quark momenta are p and p' . There is 
momentum conservation in each vertex, with all momenta inwards. 



Appendix B 

In this Appendix we show that the properties (a) to (e) in Section 1 are sufficient 
to determine (6). Conditions (d) and (e) are expressed by the Ward identities (9), (13) 
and (15) . To prove imiqueness, we suppose that there was a tensor W^j^ which could be 
added to (6) so as still to satisfy the conditions. Since the Ward identities are linear, W 
would have to satisfy them but with zero on the righthand sides. We will show this implies 
W = {). 

The most general form for W allowed by conditions (a) to (c) is 

= + Sry^r?^ + r7a/3(C^^Q" + C'-Q^) 

+QaQ(3{J''Q" + J'^Q"") + Q^Q'^Q'^Lp 

+MQ''Q''Qa,Qp + {a^(3). (SI) 

Here the tensors, C, F^"^, etc, not written out explicitly are constructed from p\^ rx (not 
using Q\ or ry^p) and scalar coefficients. The tensors constructed from Q and r] are shown 
explicitly. Note that tensors made from a single Q and no r] are not allowed by properties 
(a) to (c). A and B arc constants. 

Now use (9) (with zero on the righthand side) and (13) and equate to zero the coeffi- 
cients of the tensors which contain one 77. This gives (in the notation of (7)) 

= -AP^", = -ARi^, Da = -2AKa, T = AP.R, N = AK^, 

= -BK^, E'" = -BK", S = BP.K, S' = BR.K, = -BPp, = -BRp. {B2) 
Similarly, taking the terms with just one Q and no rj gives 
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Then taking the terms with two Q and then with three Q gives also 

J^' = J"" ^ La = M = 0. (S4) 

Finally we use the Weyl identity (15). This implies 

4:A + B + Q.D = 0, E^' + V'^ + AC'^ = 0, E'" + V"^ + = 0, AT + S + S' = 0. (S5) 

The last equation in (B5) is sufficient, with the use of (B2), to show that A — B — and 
hence that W = 0. 

FTB and JF thank CNPq and FAPESP (Brasil) for support. JCT thanks the Physics 
Department, University of Tasmania for hospitality. 
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Figure caption 

Fig.l The Feynman graphs contributing to equations (6) and (8). Doubled lines denote 
quarks, solid lines denote gluons and broken lines denote ghosts. The Feynman rules are 
in Appendix A. 



13 



